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Abstract. Neurofibromatosis type 1 (NF1) is a common auto- 
somal dominant genetic disorder caused by mutations in the 
NF1 gene. One of the hallmarks of NF1 is the high mutation 
rate in this gene. In this study, we present 127 different NF1 
mutations and 54 novel mutations detected at both the genomic 
DNA and mRNA level using a retrospective case series review. 
We found that 25.2% of these different mutations induced 
aberrant splicing. Of note, 40.6% of these splicing errors were 
caused by exonic variants. In addition, one mutation produced 
mosaicism in the post-transcriptional profile. However, studies 
investigating these splicing aberrations are limited. In order to 
better understand the pathogenicity of NF1 and to provide a 
more accurate interpretation in molecular diagnostic testing, 
combined computational analyses were employed to elucidate 
the underlying mechanisms of the variants modulating NF1 
gene splicing. 

Introduction 

Neurofibromatosis type 1 (NF1) (OMIM 162200) is a progres- 
sive autosomal dominant inherited disease and is one of the 
most widespread genetic disorders worldwide with a prevalence 
of 1 in 2500- to -3000 live births (1). The clinical characteristics 
in the NF1 diagnostic criteria include cafe-au-lait spots, neurofi- 
bromas, Lisch nodules, intertriginous freckling, typical osseous 
lesions and optic pathway gliomas (2). At least 78% of patients 
who fulfill the NIH diagnostic criteria for NF1 have NF1 gene 
mutations (3). Moreover, 5-10% of the cases are caused by a 
deletion in the NF1 gene (4). However, the positive rates of the 
NF1 mutation findings in clinical diagnostic laboratories vary 
considerably according to the proportions of the samples from 
clinically definite or suspected patients . NF1 is located on 17ql 1 .2 
and spans 28,2751 bp in length. This gene contains 60 exons and 
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encodes neurofibromin, a key component in the RAS-MAPK 
signaling pathway. The RAS-MAPK pathway regulates 
the proliferation and differentiation of neuronal cells and 
myocytes (5). Neurofibromin functions as an inhibitor of RAS 
activation and as a tumor suppressor with a central region that 
is homologous to RAS-GTPase activation proteins (GAPs) (6). 
Mutations in the NF1 gene cause a loss in neurofibromin 
function, resulting in downstream cell growth activation (7-9). 
Previous studies have reported over 1,400 different mutations 
due to the high mutation rate of the NF1 gene. A high number 
of these mutations arise as novel mutations; however, there is no 
hot spot for the pathogenic variations of NF1 (3,10,11). 

In this study, we performed a retrospective review of 
378 cases and compiled the mutations identified at both the 
genomic and mRNA level. We present 127 different mutations 
of the NF1 gene; 54 of which are novel mutations. In addition, 
the deletion of the NF1 gene was detected in 5 cases using 
fluorescence in situ hybridization (FISH) or the comparative 
genomic hybridization (CGH) array method. With the advent 
of the genomic DNA and cDNA sequencing approach, splicing 
abnormalities caused by exonic variants were captured and 
presented in our data. In addition, 7 of these 13 exonic muta- 
tions were novel. Of note, one of these mutations, c.3362A>G, 
produced mosaicism of a point mutation and mutant exon skip- 
ping at the mRNA level. 

Accurate splicing of pre-mRNA is not only controlled by 
the 5'/3' splice sites (ss), but also by other cis-acting elements, as 
well as trans-acting factors, i.e., SR proteins and heterogeneous 
nuclear ribonucleoproteins (hnRNPs). These cis-acting elements 
generally include the splicing enhancers related to exon-inclu- 
sion enhancement, splicing silencers related to exon-inclusion 
inhibition, the intronic branch point and the polypyrimidine 
tract (12,13). Although the NF1 mutation spectrum continues 
to expand, studies investigating these splicing aberrations are 
limited (14-18). Thus, integrated analyses using the bioinfor- 
matics tools were further applied to provide insight into the 
mechanisms of these splicing defects caused by exonic variants, 
as well as other intronic variants at non-consensus splice sites. 

Patients and methods 

Patients. A total of 378 cases were referred for NF1 gene testing 
in our laboratory from January, 2006 to May, 2013 and were 
recruited in this study. The subjects consisted of 338 unrelated 
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probands with clinically definite or suspected NF1 diagnosis 
and 40 family members. Consent forms were signed by the 
patients or authorized representatives. All cases underwent a 
NF1 gene-sequencing test developed in our laboratory, which 
was approved by the Ethics Committees at the University of 
Oklahoma Health Sciences Center, Oklahoma City, OK, USA. 

Mutation screening by Sanger sequencing. Genomic DNA 
was isolated from peripheral blood samples of the patients 
using the QIAamp DNA Mini kit (Qiagen, Valencia, CA, 
USA). mRNA was isolated from the peripheral blood samples 
using the QIAamp RNA Blood Mini kit (Qiagen). First-strand 
cDNA was reverse-transcribed using the Superscript III 
Reverse Transcriptase kit and random primers (both from 
Invitrogen, Carlsbad, CA, USA). PCR was performed using 
specific primers targeting the mRNA coding region of the 
NF1 gene. For confirmation, exon-specific genomic DNA 
sequencing was also performed using specific primers. Primer 
information will be provided upon request. Sanger sequencing 
was performed using the BigDye Terminator v3.1 Cycle 
Sequencing kit (Life Technologies, Foster City, CA, USA) 
and an ABI 3130x1 genetic analyzer (Life Technologies). 
Sequences were analyzed using Mutation Surveyor software 
(SoftGenetics, State College, PA, USA). 

In silico analysis. Splice Site Prediction by Neural Network 
(SSPNN; www.fruitfly.org/seq_tools/splice.html) and the 
Human Splicing Finder (HSF; www.umd.be/HSF/) were used 
to investigate the mechanisms of the splicing abnormalities 
caused by the mutations at the non-consensus splice sites. HSF 
contains its own programs and other prediction platforms, 
including the exonic splicing enhancer (ESE) finder (http:// 
rulai.cshl.edu/cgi-bin/tools/ESE3/esefinder.cgi?process=home), 
RESCUE-ESE (http://genes.mit.edu/burgelab/rescue-ese), 
FAS-ESS (http://genes.mit.edu/fas-ess), Putative Exonic 
Splicing Enhancers/Silencers (PESX) designed by Zhang and 
Chasin (30) and splicing silencer motifs designed in the study 
by Sironi et al (19). This assessment system contains 2 sets 
(HSF and SSPNN) to examine the potential splice sites, where 
1 set (HSF) was used for the potential branch points, 4 sets 
(PESx, RESCUE-ESE, ESE finder and HSF) for the ESE and 
3 sets (PESx, FAS-ESS and splicing silencer motifs) for the 
exonic splicing silencer (ESS). The query sequences were 
obtained from the normal and mutated sequences. The SSPNN, 
HSF and ESE finder provided the scores to value the strength 
of the splicing-relative sequence motifs using corresponding 
weight matrices. PolyPhen2 (http://genetics.bwh.harvard.edu/ 
pph2) and SIFT (http://sift.bii.a-star.edu.sg) were applied to 
predict the potential effect of an amino acid substitution on the 
structure and function of the NF1 protein. BLASTP was used 
to align the NF1 protein sequences from the multiple species, 
including human, chimpanzee, gorilla, cat, dog, mouse, rat, 
cattle, chicken and zebrafish. 

Real-time PCR. The alternative post-transcriptional profiles in 
one specific case were produced using the SYBR Master Mix 
(Life Technologies) and ABI PRISM 7000 Sequencing Detection 
System (Life Technologies). The ACTB gene encoding (3-actin 
was used for normalization. For ACTB, the forward primer was 
5-AGCTCCTCCCTGGAGAAGAG-3' and the reverse primer 



was 5'-AGCACTGTGTTGGCGTACA-3'. For NF1, the forward 
primer was 5'-GATGTAAAATGTCTTACAAG-3' and the 
reverse primer was 5'-CTGCCACCTGTTTGCGCACT-3'. 
Amplicons targeting on the NF1 and ACTB genes were confirmed 
using Sanger sequencing. Real-time PCR was performed in 
triplicate using 5, 2.5, 1.25 and 0.625 ng cDNA. 

Results 

Mutation spectrum. Mutational screening of the NF1 gene 
was performed on samples obtained from 378 clinically diag- 
nosed or individuals suspected of having NF1 . The mutation 
nomenclature was based on the NCBI reference NM_000267.3. 
The exon number was given according to the conventional 
rule used in the NF1 testing community and previous litera- 
tures (10,11,17,18). The mutations were confirmed using the 
Biobase (HGMD professional version database) and Leiden 
Open Variation Database (LOVD) to determine the recur- 
rence. In addition, the missense mutations were searched in 
1000 Genomes project, dbSNP and Exome Variant Server to 
rule out normal variants. NF1 mutations were identified in 
169 out of 378 cases; 127 different mutations were observed 
(Tables I and II). Of these mutations, 54 mutations were novel, 
of which, 23 were frameshift mutations, 10 were splicing 
defects, 7 were nonsense mutations and 14 were missense muta- 
tions (Table III). The mutations affected almost all exons apart 
from exon 4c, 14, 23.1, 35, 38 and 49 in the mutation spectrum 
of the patients, which was consistent with the finding of no hot 
spot mutations in previous studies (3,10,1 1). The nonsense muta- 
tions were the most common molecular defects found in this 
study (33/127), followed by the splice-site mutations (32/127) 
and missense mutations (27/127). In the group of frameshift 
mutations, deletion (n=23) was prone to occurring compared 
to insertion/duplication (n=8) and indels (n=4), but most of the 
insertion/duplication and all of the indels were novel mutations. 

Only 50% of the splicing defects disrupted the conserved 
GT/AG or AT/AC dinucleotides of the splice sites in this study 
(Table II). By contrast, 3 intronic mutations at non-consensus 
splice sites (Table IV, subgroup I) generated cryptic 5'ss or 
3'ss, resulting in the insertion into the mRNA. Thirteen exonic 
variants, that were 40.6% of the splicing defects, induced exon 
skipping or aberrant exons instead of a point mutation based 
on the genomic DNA and cDNA sequencing methods. Among 
these exonic mutations, c.1466A>G and c.1885G>A (Table IV, 
subgroup II) induced aberrant exons with a deletion by gener- 
ating cryptic 5'ss or 3'ss, respectively. In particular, the 1 185G>A 
mutation was a silent mutation at the genomic DNA level. The 
other 4 exonic sequence alterations, which resulted in exon 
skipping, were identified as substitutions at the last nucleotide 
position of exon 8, 20 and 29, and the last second nucleotide of 
exon 11 (Table IV, subgroup III). The remaining 7 exonic muta- 
tions, which caused exon skipping, did not perturb the natural 
3'/5'ss or create cryptic splice sites (Table IV, subgroup IV). 

Of note, one exonic mutation, c.3362A>G, produced 
mosaicism of E1121G and exon 20 skipping at the mRNA 
level (Fig. 1). This germline mutation in the proband was inher- 
ited from his father, and the post-transcriptional mosaicism was 
presented in the mRNA of both patients. However, the missense 
mutation of the son showed a lower signal intensity using Sanger 
sequencing (Fig. 2). The mosaicism of the post-transcriptional 
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Table I. The 127-mutation spectrum apart from the 


splicing 


abnormalities. 






cDNA 


Protein 


Exon 


cDNA 


Protein 


Exon 


Missense mutations 












1A>G 


MetlVal 


1 


3827G>A 


Argl276Gln 


22 


410C>G 


Serl37Cys 


4a 


4172G>C 


Argl391Thr 


24 


1241T>G 


Leu414Arg 


9 


4288A>T 


Asnl430Tyr 


25 


1646T>C 


Leu549Pro 


11 


4306A>G 


Lysl436Glu 


25 


2350 T>C 


Trp784Arg 


15 


4306A>C 


Lysl436Gln 


25 


2975T>A 


Met992Lys 


17 


4318A>G 


Metl440Val 


25 


3104T>C 


Metl035Thr 


18 


4493G>A 


Glyl498Glu 


26 


3142T>G 


Trpl048Gly 


19a 


5131G>C 


Alal711Pro 


28 




Ala 1071 Vrf\ 


19b 


5425C>T 


Argl809Cys 


29 


3295 A>G 


Lysl099Glu 


19b 


5489C>T 


Prol830Leu 


29 


3296 A>C 


Lysl099Thr 


19b 


5498T>C 


Leul833Pro 


29 


3362A>G 


Glull21Gly 


20 


5855C>T 


Alal952Val 


31 


3494T>C 


Ilell65Thr 


20 


7106T>C 


Leu2369Ser 


39 


3752A>C 


Hisl251Pro 


11 








Nonsense mutations 












569T>G 


Leul90* 


4b 


4084C>T 


Argl362* 


23.2 


574C>T 


Argl92* 


4b 


4107C>G 


Tyrl369* 


23.2 


586G>T 


Glul96* 


4b 


4267A>T 


Lysl423* 


24 


668G>A 


Trp223* 


5 


4537C>T 


Argl513* 


27a 


1238C>G 


Ser413* 


9 


5264C>G 


Serl755* 


29 


1246C>T 


Arg416* 


9 


5401C>T 


Glnl80r 


29 


1275G>A 


Trp425* 


10a 


5708T>G 


Leul903* 


30 


1318C>T 


Arg440* 


10a 


5839C>T 


Argl947* 


31 


1381C>T 


Arg461* 


10a 


5941C>T 


Glnl981* 


31 


1754T>A 


Leu585* 


12a 


6243C>G 


Tyr2081* 


33 


2041C>T 


Arg681* 


13 


6709C>T 


Arg2237 * 


36 


Zj jZu>A 


irp / 54 


15 


7285C>T 


Arg2429* 


41 


2446C>T 


Arg816* 


16 


7486C>T 


Arg2496* 


42 


3049C>T 


Glnl017* 


18 


7843C>T 


Gln2615* 


45 


3826C>T 


Argl276* 


22 


7993C>T 


Gln2665* 


46 


4006C>T 


Glnl336* 


15.1 


8072G>A 


'!■ 1 / 111 * 

lrp2691 


47 


4066G>T 


Glul356* 


132 








Frameshift mutations 












118delA 


Lys40Argfs*4 


1 


4498_4505delTATCTTTC 


Tyrl500_Ilel501delfs*6 


26 


154delT 


Ser52euLfs*4 


1 


4688_4689insAA 


Phel563Leufs*5 


27b 


499_502delTGTT 


Cysl67Glnfs*12 


4b 


4810dupT 


Tyrl604Leufs*16 


28 


802_803delCCinsG 


Pro268Aspfs*13 


6 


4906dupA 


Aspl636Argfs*5 


28 


1010_1021dell2 


Glu337_Ser340del 


7 


4930_4937indels 


Vall644Serfs*3 


28 


1541_1542delAG 


Gln514Argfs*43 


10c 


5010delG 


Lysl670Asnfs*7 


28 


1664_1667delTAGA 


556Aspdelfs*13 


11 


5040_5043delAGGTinsTA 


Lysl680Asnfs*16 


28 


1756_1759delACTA 


Phe586delfs*19 


12a 


5909dupC 


Ilel971Tyrfs*17 


31 


1882delT 


Tyr628Thrfs*3 


12b 


6791dupA 


Tyr2264* 


37 


1908delT 


Ser636Valfs*52 


12b 


7096_7101delAACTTT 


Asn2366_Phe2367del 


39 


2032_2033delCCinsA 


Pro678Lysfs*10 


13 


7125delA 


Tyr2375Thrfs*20 


39 


2342dupA 


His781Glnfs*13 


15 


7160delG 


Arg2387Lysfs*10 


40 


2984_2988delTGGTC 


Leu995Glnfs 24 


17 


7267dupA 


Thr2423Asnfs*4 


41 


3054delT 


Leul018* 


1 o 


7^88 7^84-1-1 "?Hp110 


\J lUSJlU W 11 


zL1 

T- 1 


3108_3109ins23 


Lysl036Thrfs*8 


18 


7510delG 


Asp2504Thrfs*23 


42 


3852_3854delAAT 


Llel284fs 


22 


7581_7582delAT 


Ser2528Glnfs*21 


43 


4312_4314delGAA 


Glul438del 


25 


7726delT 


Ser2576Glnfs*27 


44 


4418_4419delAT 


Hisl473Glnfs*7 


26 









Novel mutations were denoted in bold text. 
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Table II. The mutation spectrum of the splicing abnormalities. 



Splicing mutations at the consensus splice sites 



Mutation 


iva 


cDlNA errect 




Effect on splice site 


60+ldelG 


IVSl+ldelG 


Unknown 




Inactive 5'ss 


204+lG>A 


IVS2+1G>A 


100_204dell05 




Cryptic 5'ss 


205-2A>C 


IVS2-2A>C 


AE3 




Inactive 3'ss 


889-2A>G 


IVS6-2A>G 


AE7 




Inactive 3'ss 


889-lG>A 


IVS6-1G>A 


AE7 




Inactive 3'ss 


1642-1G>A 


IVS10c-lG>A 


AE11 




Inactive 3'ss 


2409+ 1G> A 


IVS15+1G>A 


AE15 




Inactive 5'ss 


2850+ 1G> A 


IVS16+1G>A 


AE16 




Inactive 5'ss 


3114-2A>G 


IVS18-2A>G 


AE19a 




Inactive 3'ss 


3708+2T>A 


IVS21+2T>A 


AE21 




Inactive 5'ss 


3709-2A>G 


IVS21-2A>G 


3709_37 1 8delGATGAACTAG 




Cryptic 3'ss 


4367+ 1G>A 


IVS25+1G>A 


AE25 




Inactive 5'ss 


6579+ 1G>A 


IVS34+1G>A 


AE34 




Inactive 5'ss 


6858+lG>A 


IVS37+1G>A 


AE37 




Inactive 5'ss 


7676-2A>G 


IVS43-2A>G 


AE44 




Inactive 3'ss 


8098-lG>A 


IVS47-1G>A 


AE48 




Inactive 3'ss 


Intronic mutations at non-consensus splice sites 


Mutation 


Related IVS 


cDNA effect 




Cryptic splice site 


1260+3A>T 


IVS9+3A>T 


1261_1262insGTTAGTCCAAAAG 




Cryptic 5'ss 


2410-18C>G 


IVS15-18C>G 


2410_2411insl7 




Cryptic 3'ss 


5944-5A>G 


IVS31-5A>G 


5943_5944insCTAG 




Cryptic 3'ss 


Exonic mutations 


Mutation 


Related exon 


cDNA effect 




Cryptic splice site 


910OT 


In E7 


AE7 




No 


972T>A 


In E7 


AE7 




No 


1185G>A 


Last NT of E8 


AE8 




No 


1466A>G 


In ElOb 


1466_1572del62 




Cryptic 5'ss 


1720 A>G 


Last second NT of El 1 


AE11 




No 


1885G>A 


In E12b 


1846_1886del 




Cryptic 3'ss 


3362A>G 


In E20 


AE20 




No 


3467A>G 


In E20 


AE20 




No 


3496G>A 


Last NT of E20 


AE20 




No 


5546G>A 


Last NT of E29 


AE29 




No 


6792C>A 


In E37 


AE37 




No 


6792C>G 


In E37 


AE37 




No 


7694delC 


InE44 


AE44 




No 


Novel mutations were denoted in bold text, ss, splice sites. 


Table III. The number of each type of abnormality in 


127 mutations and the frequency of the novel mutations. 




Missense Nonsense 


Deletion Ins/Dup 


Indels 


Splicing defect 


Total no. 


27 33 


23 8 


4 


32 


Novel mutations 


14 7 


13 6 


4 


10 
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CTT S G A G A T G AAA 



Proband's genomic DNA c. 3362A>G 




T CT TCAACT T CA 



Proband's cDNA c. 3362A>G (on reverse strand) 




Figure 1 . The Sanger sequencing data showed a germline mutation (c .3362 A>G) 
resulting in the missense mutation, E1121G, and exon 20 skipping coexisting 
at the mRNA level. 



AT C T TCAACT T C 

Father's cDNA | Mutation on reverse strand 




Figure 2. Comparison of the signal strength of the missense mutation in the 
mosaic cases. 



profile was further examined using real-time PCR with cDNA 
obtained from this family (Fig. 3). The sample from the mother 
was used as the wild-type sample in this assay. The skipping 
rate of the mutant exon in the samplefrom the son was higher 
compared to the rate in the sample from the father. 

In silico analysis. To better understand the underlying mecha- 
nisms of these 16 unusual splicing errors, 7 computational 
tools were employed to examine how these mutations affected 
the splicing-relative sequence motifs. The comparison was 
made between the results of the prediction on the normal and 
mutated sequences (Table IV). SSPNN predicted a deletion of 
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Figure 3. Real-time PCR profiles of exon 20 skipping in the family with the 
germline mutation c.3362A>G found in the son and father. As the template 
amount decreased, the cycle number increased in the amplification plot. The 
data indicated that exon 20 template was less in the sample of the son. 



the authentic acceptor site (3'ss) caused by 2410-18C>G, and a 
decrease in the strength of the authentic splice sites caused by 
1260+3A>G and 5944-5A>G in the subgroup I. By contrast, 
the cryptic splice sites generated by these 3 mutations and 
the subgroup II mutations were given high scores by HSF or 
SSPNN. Substitutions in the subgroup III all abolished the 
authentic donor sites (5'ss) in the SSPNN evaluation, compared 
to the prediction of the strength reduction by HSF. No cryptic 
splice sites were predicted in these substitution sequences and 
the subgroup IV sequences. 

Further investigation on other splicing-relative sequence 
motifs revealed that some cis-acting factors were altered in 
these mutated sequences. In the subgroup I sequences, the 
branch-point motif was deleted and a putative ESE was gener- 
ated by 2410-18C>G; a putative ESS or an abolished ESS was 
predicted to accompany the strength-decreased donor site or 
acceptor site in the mutated sequences of 1260+3A>G and 
5944-5A>G, respectively. In the subgroup II sequences, the ESE 
was deleted or the strength of the ESE was decreased simulta- 
neously with the strong new splice site generated by 1466A>G 
and 1885G>A, respectively. The original ESE was deleted or no 
ESE was embedded in the subgroup III sequences. The archi- 
tectural alterations of the splicing-regulatory elements were 
more complex in the subgroup IV sequences. In general, the 
decreased strength of the ESE or decreased ratio of the ESE/ 
ESS was presented. An abolished ESS was also predicted in 
mutated exon 7 and 37, but no abolished ESE was predicted. 

Discussion 

NFl is a multisystem genetic disorder with extreme diversity 
of clinical expression (3,10,11,20,21,22). A clear correlation 
of genotype and phenotype has been previously demonstrated 
in only two types of mutations (23). Patients with an NFl 
microdeletion have more severe clinical characteristics (4,23). 
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We found 5 cases with an NF1 gene deletion in patients 
with various severe conditions, such as developmental delay, 
seizures, and early onset skin/subcutaneous tissue disorders. 
One 4 year-old patient harboring a novel splice site mutation 
(60+ldelG) also exhibited the developmental delay. This alter- 
nation, which induced skipping in exon 1 may cause the same 
effect as an NF1 deletion at the protein level. Another mutation, 
c.2970-2972delAAT, in exon 17, which was not shown in our 
study, has been reported to be associated with the absence of 
cutaneous neurofibromas (11,24). However, due to the lack of 
detailed clinical information, our investigation of the correla- 
tion between genotype and phenotype was limited. The present 
study focused largely on the mutation spectrum, particularly 
the splicing errors. 

The mutation rate of the NF1 gene is one of the highest 
reported in the human genome. We presented a mutation spec- 
trum of the NF1 gene with 127 different variants in this study. 
In the summarized data, exon 19b, 25 and 29 harbored more 
missense mutations, exon 4b, 10a and 23.2 appeared to have 
more nonsense mutations, and the highest rate of frameshift 
mutations was found in exon 28. In addition, exon 7, 20 and 37 
were prone to splicing errors. A total of 54 of the 129 different 
mutations were novel mutations that were categorized as either 
frameshift mutations, splicing defects, nonsense mutations or 
missense mutations. With the exception of missense mutations, 
the other 3 types of mutations are considered highly likely to be 
deleterious. Patients with NF1 with a missense mutation have 
a lower incidence of multiple neurofibromas and plexiform 
neurofibromas compared to patients with a different type muta- 
tion; In addition, it is also true that no evidently milder NF1 
phenotype was concluded to be distinctly associated with a 
missense mutation (23). The in silico analysis of the functional 
consequence of these novel missense mutations was performed 
using the Polyphen2 and SIFT program. A total of 12 of these 
14 novel missense mutations were predicted to be potentially 
damaging. Although P1830L and A1952V were determined to 
be benign, these mutations and other novel missense mutations 
altered amino acids that were conserved across the different 
species according to the results obtained from an orthologous 
alignment using BLASTP. It is widely believed that mutations 
in a highly conserved region may result in a functional change. 
Moreover, we found 2 second missense variants, M645V and 
I1658V, with nonsense mutations. These 2 variants can be 
treated as neutral polymorphism after the parental study. 

In the present study, 25.2% of the different mutations induced 
aberrant splicing and 50% of these splicing errors were caused 
by exonic mutations and intronic mutations at non-consensus 
splice sites. These results are consistent with those of previous 
findings, showing that the NF1 gene is susceptible to having 
splicing errors in post-transcription (14,15). Accurate splice 
site recognition is critical in pre-mRNA splicing (25). This 
process is a coordinated program involving the strong splice 
sites, correct splicing regulatory elements (SRE) embedded in 
the genome and the associated proteins. Bioinformatics assess- 
ments focused on revealing how the authentic splice sites and 
the hidden sequence motifs of these SREs were interfered by 
these exonic and intronic variants at non-consensus splice sites 
in this study. We found that the change in these splice sites acted 
by tethering the alteration of the ESE, ESS and other cis-acting 
elements, which resulted in aberrant splicing. The error-prone 



splicing occurred in the subgroups I and II in which the 
high-score cryptic 5'/3'ss or the next AG/GT with a higher score, 
as the case of 1260+3A>G, replaced the strength-decreased or 
lower-scored authentic ones in the new microenvironment of 
splicing regulatory elements. For example, the 2410-18C>G 
mutation generated a higher score cryptic splice site while 
forming a putative ESE and abolishing the original branch 
point, and then 17 base pairs were inserted into the mRNA 
as the consequences of a strong cryptic splice site in coor- 
dination with the gain and/or loss of SREs. Importantly, the 
2410-16A>G, 2410-15A>G, 2410-12T>G mutations have been 
previously documented in the Biobase HGMD database. Given 
2410-18C>G, it is an obvious sign to alert that this cluster near 
the intron 15/exon 16 boundary is a splicing-aberration harbor. 

In some cases, the disruption of ESE elements was the 
principal cause of the splicing error due to the consequences 
of the reduced splicing enhancement activity (16,26,27). As 
is known, not all of the substitutions at the last nucleotide 
position of an exon will cause exon skipping, although the 
mutants are also predicted to delete the authentic splice site or 
decrease its strength. The loss of ESE motifs may explain why 
the exon skipping was caused by the exonic mutations in the 
subgroup III. An exception in this case was exon 29 that had 
no ESE motif. Exon 29 skipping resulted from the weakened 
donor site caused by the mutation 5546G>A and the lack of the 
ESE to support the splice-site recognition. 

There was no cryptic splice site generated and no strength 
reduction of the authentic splice sites in the subgroup IV 
sequences. However, the acceptor site strength was character- 
ized to be an important and sensitive parameter in splice-site 
recognition (28). The low-score exons, such as exon 7 and 37 were 
observed to have more splicing defects in this study. In addition, 
our in silico analysis also showed more complex changes in the 
ESE and ESS. Not only was the ESE demonstrating a decrease 
in strength, but it also showed that weaker ESS, decreased 
ratio of ESE/ESS and abolished ESS were predicted to be the 
architectural weakness in the exon definition, resulting in exon 
exclusion in the subgroup IV sequences. The ratio of ESE/ESS 
was important for exon recognition and intron identification in 
the complexity of splicing. Several studies have observed that 
a higher density of ESEs was in the exons compared to the 
introns and vice versa for ESSs (12,29,30). The ratio decreasing 
of ESE/ESS will break the delicate balance of SREs, and will 
be prone to exon exclusions. Although enhancers and silencers 
have apparently opposite effects as suggested by their terms, the 
Composite Exonic Regulatory Elements of Splicing (CERES) 
has already been proposed when accumulating evidence has 
suggested ESE and ESS shared additional properties (31,32). 
The findings of the weaker or abolished ESS caused by the 
mutations in our analysis reflected the overlapping function 
of these 2 elements. The SREs became more critical in deter- 
mining the splice-site recognition in these cases. 

When the trans-acting factors navigated the new landscape 
in which the mutation plays a make-or-break role, the conse- 
quences of the competition and coordination in the splicing 
process is more evident in the case with the 3362A>G muta- 
tion, in which the germline mutation decreased the ratio of 
the ESE/ESS. The splicing in the mutated microenvironment 
resulted in the missense mutation and exon skipping coex- 
isting in the mRNA of an 11-year-old boy and his father with 
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different rates. Furthermore, real-time PCR clearly confirmed 
this mosaicism. Both patients have multiple cafe-au-lait spots 
but no other profoundly different and NFl-related clinical 
features. The higher skipping rate found in the sample of the 
son may be caused by individual's genetic variability. 

Taken together, this study presents 54 novel NF1 mutations 
and reveals the high frequency of the unusual splicing defects in 
the pathogenicity of NF1. Integrated analyses using the bioin- 
formatics tools provided insight in order to better understand the 
underlying mechanisms of these splicing errors. In particular, 
as the NF1 gene was susceptible to having aberrant splicing 
caused by exonic variants, including the silent mutation, such 
as 1 185G>A, this result underscored the large consequences 
of NF1 gene testing at both the genomic and mRNA levels. 
In addition, the mutation data may contribute information to 
the ongoing antisense therapeutics for NF1 caused by intronic 
mutations (33), thus shedding light on a targeted treatment to 
restore NF1 gene function. 
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